THE HEUN EQUATION AND THE 
CALOGERO-MOSER-SUTHERLAND SYSTEM V: GENERALIZED 
DARBOUX TRANSFORMATIONS 
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Abstract. We obtain isomonodromic transformations for Heun's equation by gen- 
eralizing Darboux transformation, and we find pairs and triplets of Heun's equation 
which have the same monodromy structure. By composing generalized Darboux 
transformations, we establish a new construction of the commuting operator which 
ensures finite-gap property. As an application, we prove conjectures in part III. 



1. Introduction 

It was shown in [Hj that some pairs of Schrodinger operators are isomonodromic. 
Set 

(1.1) H, = -^ + 6p{x), 

(1.2) H2 = + 2p{x) + 2p{x + cui) + 2p{x + UJ2), 

where p{x) is the Weierstrass p-function with periods {2ui, 2^3) and a'2 = —ivi — lu^. 
Now we consider eigenfunctions of Hi (resp. H2) with the eigenvalue E. Set 

Ei{x, E) =Qp{xf + 'iEp{x) + E^- 9(72/4, 

H2(x, E) ={E - 3e3)p(x) + {E - 3e2)p(x + Ui) 

+ {E- 3ei)p(x + UJ2) + E^~ 3^2/2, 

Q{E) ={E^ - 3g2)iE - ?>e{){E - 3e2)(^ - 363), 



r^/^Q{E)dx 

Ak{x,E) =^/^k{x,E)exp / — -— , 

J ^k{x,E) 



(A; = 1,2), 



where Cj = p{u!i) {i = 1,2,3) and g2 = — 4(eie2 + 6263 -|- 6361). Then it was shown 
that the functions Ai(a;, E) and Ai(— x, E) (resp. A2{x, E) and A2(— x, E)) are eigen- 
functions of Hi (resp. H2) with the eigenvalue E, and they satisfy 

A,(±(x + 2u;,),i^) = A.(x,E)exp I t^- I Z^Elh±^£^M^dE 




-Q{E) 

for k = 1,2 and z = 1, 2, 3, where r]i = ({uJi) and ({x) is the Weierstrass zeta function. 
Hence the monodromy of eigenfunctions of Hi with the eigenvalue E coincides with 
that of eigenfunctions of H2. 

In this paper, we investigate this phenomena by Darboux transformation and gen- 
eralized Darboux transformation. Let (f)o{x) be an eigenfunction of the operator 
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H = —cP/dx"^ + q{x) with an eigenvalue Eq, i.e. 



dx"^ 



For this case, the potential is written as = {(p'^i^x) / (f)Q{x)y + {(pQ^x) / (t>o{x)Y + 
Eq. Set L = d/dx - (p'^ix) / Mx) and H = -d^dx^ + q{x) - 2(0(j(x)/0o(x))'. Then 
we have 

HL = LH. 

Hence, if is an eigenfunction of the operator H with the eigenvalue E, then L(f){x) 
is an eigenfunction of the operator H with the eigenvalue E. This transformation is 
called the Darboux transformation. We generalize the operator L to be the differential 
operator of higher order, and we call it the generalized Darboux transformation. 

The Schrodinger operator we consider in this paper is the Hamiltonian of the BCi 
Inozemtsev model, which is written as 

(1.3) H^ioA,hM ^ + J2 kih + l)pix + cu,), 

where uq = 0. The potential of this operator is called the Treibich-Verdier potential, 
because Treibich and Verdier |13j found and showed that, if li G Z>o for all i G 
{0, 1, 2, 3}, then it is an algebro-geometric finite-gap potential. For further results on 
this subject, see [21 El 112] • The algebro-geometric finite-gap property cause the 
possibility for calculation of eigenfunction and monodromy of the operator if ('o>'i>'2,'3)_ 
Let f{x) be an eigenfunction of the operator if ('o>'i>'2,'3) -^ith the eigenvalue E, 
namely, 

(1.4) (^--^ + J2 hik + l)p{x + ^.) j fix) = Ef{x). 

Then this equation is an elliptic representation of Heun's equation. Here Heun's 
equation is the standard canonical form of a Fuchsian equation with four singularities 
(see jl]). Thus, solving Heun's equation is equivalent to studying eigenvalues and 
eigenfunctions of the Hamiltonian of the BCi Inozemtsev model. 

We now describe the main result of this paper. Let ctj be a number such that 
a-i = —li or ttj = /j + 1 for each i G {0, 1,2,3}. Set d = — Yli=o'^i/'^ assume 
d G Z>o. Then there exists a differential operator L of order d + 1 which satisfies 

jj{ao+d,ai+d,a2+d,a3+d) _ j^j^{lo,h,h,h) 

Note that Khare and Sukhatme [3^ essentially established this result for the case 
d = 0, that is the case of original Darboux transformation. It follows immediately 
that, if 0(a;) is an eigenfunction of the operator H'^^ohMM -^^ith an eigenvalue E^ then 
L(j)[x) is an eigenfunction of the operator ii("o+d,«i+Q!,a2+(i,a3+d) -^^ith the eigenvalue 
E. Since all coefficients of the operator L with respect to the differential {d/dxY 
{k = 0, . . . ,d + 1) is shown to be doubly-periodic, the operator L preserves the 
data of monodromy. Hence the operators H^^o,h,h,i3) g^^^ j^{ao+d,ai+d,a2+d,a3+d) 

isomonodromic, and isospectral, because boundary condition for spectral problem 
is characterized by monodromy. Note that the condition d G Z>o corresponds to 
quasi-solvability of the operator ii('0''i'^2,;3)_ 
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For the case that Iq, h, I2, h are all integers, there exists an operator //('o-'i'^a-'s) g^ch 
that the pair ijCo-'i'^a-'a) jyihAMM jg connected by isomonodromic transforma- 
tion. In some cases, they are self-dual. For example, the operator Hi(= ff(2.o,o,o)-j 
Eq. ljl.lj) is connected to the operator H2{= H^^'^'^'^^) in Eq. (jl.2j) by the transforma- 
tion L = d/dx — p'{x)/{2{p{x) — ei)) — p' {x) / {2{p{x) — 62)), i.e. we have 

^(1,1,1,0)^ ^ (2,0,0,0) _ 

For the case that /q + 1/2, /i + 1/2, Z2 + 1/2, /s + 1/2 are all integers, there exists 

,,(1) ,(1) ,(1) ,(1), ,j(2) ,(2) ,(2) ,(2), I] 1 1 1 \ 

two operators if^'o ''1 ''2 ''3 ) and H^'-o ''1 ''2 .'3 ) guch that the triplet _f/'i'0''i''2,(3j^ 
^('0 ''1 ''2 ''3 ) and if'-'o ''1 ''2 ''3 ) ig connected by isomonodromic transformations. 

In the paper finite-gap property of the operator H^^o,h,h,h) ^^j^g ^ase /q, h, I2, h ^ 
Z>o is studied (see also [121 US 12 E])- Especially, a differential operator A of odd 
order which commutes with ff('0''i''2>^3) jg constructed. In this paper, we propose a 
new method for construction of the commuting operator by composing four general- 
ized Darboux transformations. Note that each generalized Darboux transformation is 
written explicitly. To show that the commuting operator constructed by composing 
four generalized Darboux transformations coincides with the one defined in P|, we 
need a discussion which will be done in section [B] As an application, we prove con- 
jectures in [8j. Namely, we establish that the polynomial defined by quasi-solvability 
coincides with the one defined by using the doubly-periodic function S(x, E) which 
is written as a product of two eigenfunctions. We also prove that the commuting 
operator is characterized by annihilating the spaces of quasi-solvability. As for the 
isomonodromic pair _f/'('o,'i,'2,i3) and i7('o.'i>'2,;3)^ jg ghown that some functions related 
to the monodromy of if(^o>^i''2.'3) coincide with that of H^^o,h,h,h) _ 

This paper is organized as follows. In section |2l we review the connection between 
quasi-solvability and generalized Darboux transformation which was essentially done 
in pp. In section 121 we consider generalized Darboux transformations for the case 
of Heun's equation. In section |^ (resp. section |2)), we investigate isomonodromic 
transformations for the case lo,li,l2,h ^ ^ (resp. lo,li,l2,h G Z + 1/2). In section 
El we construct a differential operator of odd order which commutes with ij('0''i''2.'3) 
and investigate it. 

2. Quasi-solvability and generalized Darboux transformation 

We review the relationship between the quasi-solvability and the generalized Dar- 
boux transformation. 

We set H = —d'^/dx^ + q{x). Let be a positive integer. If the operator H 
preserve a n-dimensional space U of functions, then the operator H is called quasi- 
solvable. Then there exists a basis (/i(x), . . . , /„(x)) of the invariant space such 
that Hfj{x) = J2i^i,jfii^) f*^^ some constants (1 < i,j < n). Let PH,uif) 
be the characteristic polynomial of the operator H on the space U . Then the set 
{E\Ph,u{.E) = 0} coincides with the set of eigenvalues of of the operator H on 
the space U. Then the model is partially solved, and this is an origin of "quasi- 
solvability". For the space U, there exists a monic differential operator of order n 
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such that Lf{x) = for all f{x) G U. It is determined uniquely and written as 



f2ix) ±}2{X) 



I _d_y 

K dx) 



-1 



\ 



U{x) ±Ux) 

1 



dx 



Proposition 2.1. (^c./. [I J Assume that the operator H = —(P/dx^ + q{x) preserve 
a n-dimensional space U of functions. Let L be the differential operator written as 
Eg. \2. 1\) which annihilate the functions in U. Set H = —d'^/dx'^ + q{x) + 2ci{x). Then 
we have 

HL = LH. 

Proof. By a direct calculation, it follows that the order of the differential operator 
HL — LH is at most n — 1. Assume that HL — LH ^ 0. We denote the order by k. 
Then the dimension of solutions to the differential equation [HL — LH)f{x) = is 
k. Let g{x) G U. Since the operator H preserve the space U, we have Hg{x) G U. 
By definition of the space U, we have Lg{x) = and LHg{x) = 0. Hence we have 
{HL — LH)g{x) = for g[x) G U, but it contradicts to that the dimension of solutions 
is k{< n — 1). Therefore we obtain HL = LH. □ 

We consider the case n = 1, Let 4'o{x) be a non-zero function in U. Then U = 
C(f)o{x), and the operator which annihilate 4>o{x) is written as L = d/ dx—cp'^^x) / (p^i^x). 
The operator H is written a.s H = H — 2{(I)'q{x) / (f)Q{x))' . Hence the proposition 
reproduce the Darboux transformation. In this sense, the transformation in the 
proposition may be called the generalized Darboux transformation. 



3. Generalized Darboux transformation for Heun's equation 

In this section, we apply Proposition 12. II for Heun's equation. For this purpose, we 
recall quasi-solvability of Heun's equation. 

Proposition 3.1. 7, Proposition 5.1] Let ai he a number such that Ui = —U or 
Oii = U + 1 for each i G {0, 1, 2, 3}. Set d = — XlLo "^^/^ ^'^^ assume d G Z>o. Let 



00,01,02,03 ^6 + ^-dimensional space spanned by 
(3.1) 

Then the operator /7('0''i''2,;3) ^^gg Eq.^T^) preserves the space Va^^^ai,c 



{(p(x) - e,r'\p{x) - e,r'\p{x) - e3)"^/V(^)"}„=o,., 



Set 2 = p(x), ^{z) = (z-ei)°i/2(^_e2)"2/2(^;2_ 63)^3/2^ and H^^ohMM = ^z)'^ o 
}j{io,h,i2,h) o <l>(z). Proposition 13. II is proved by showing that the operator H^^ohMM) 
preserve the space spanned by {z — 62)^' (r = 0, . . . , (i). For details, see the proof of 
[H Proposition 5.1]. 

Now we calculate the differential operator which annihilate the space Kto,ai ,02,03- 

Proposition 3.2. The monic differential operator of order d+1 which annihilate the 
space VaQ,ai,a2,az is Written as 



(3.2) 



p'ixr+'^pix)) o 



^00,01,02,03 



1 d 
p'{x) dx 



d+1 



o$(p(a;)) 



-1 
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We write the operator La^^ai,a2,a3 o-s 

(3.3) /^ao,ai,«.,«3 - y^^j + 1^ Ci[x) y^^j 

Then 



d + l / A 2ai + d 
(3-4) ci{x) = —\2^ 



If i is even (resp. odd), then Ci{x) is expressed as Ci{x) = Ri{p{x)) (resp. Cj(x) = 
Ri{p{x))p\x)), where Ri{z) is a rational function in z. 

Proof. It is trivial that the operator {d/dzY^^ annihilate the space spanned by 
(r = 0, . . . , (i), and the operator $(z)o((i/(iz)^+^o<|)(z)~^ annihilate the space spanned 
by ^{z)z'' (r = 0, . . . , d). Write 

Then £1(2;) = —J2'i=iiid + l)ai)/{2{z — Cj)), and all coefficients Ci{z) are rational 
functions in z. By the transformation z = p{x), the monic operator 

1 d^'^' 



(3.5) L^,,a,,a,,as = p'(x)'^+^$(p(x)) O ( ) O <l>{p{x))-' 



, \ d+l d+l , , s d+l-i 

annihilate the space V^o,ai,Q2,03- Write {{l/p'{x)){d/dx)y = J2]=oh(^)(d/dxy . If j 

is even (resp. odd), then bj{x) is expressed as rj{p{x)) (resp. rj{p{x))p'{x)), where 
rj{z) is a rational function in z. Combining with the relation p'{xy = A{p{x) — 
ei)(p(a;)— e2)(p(a;) — 63), it follows that, if i is even (resp. odd), then Cj(x) is expressed 
as Cj(x) = Ri{p{x)) (resp. Ci{x) = Ri{p{x))p'{x)), where Ri{z) is a rational function 
in z. We now calculate Ci{x). By Eqs. (j3.5t IA.2|) . we have 

□ 

Theorem 3.3. Let &e a number such that ai = —U or ai = U + 1 for each 
i G {0,1,2,3}. Set d = — X]i=o'^*/2 ^'^^ assume d G Z>o. Let Lao,ai,a2,a3 be the 
operator defined in Proposition \'J/A Then we have 

(o f\\ TT{ao+d,ai+d,a2+d,a3+d) j _ j TJ [lo ,h M ,h) 

v'-'''-'/ ao ,ai ,a2 ,03 ao, 0.1,02,013 

Proof. It follows from Propositions 12. II and 13.21 that 

( TT{h,h,l2,l3) \ '^J lrp\\T_ _r fr{lo,h,l2,l3) 

V ~ ^iK"^ J J 00,01,02,03 -'-'00,01,02,03-'^-'^ ) 
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where Ci(x) is defined in Eq. ()3.4|) . By Eq. ()A.2|) we have 

= -{d + l)(2(ai + a2 + as) + Sd)p{x) + + l)(2ai + d)p{x + tOi). 

i=l 

Since = —li or = /j + 1, we have Zj(Zj + 1) = — + 1). Hence we ob- 
tain ('0''i''2''3) -\-2c[{x) = f/'(°o-l,ai-l,a2-l,a3-l) -|-2c']^(x) = ffi'^o+d,ai+d,a2+d,a3+d) 

Eq.dSSl). □ 

We consider the converse relation to Eq. (jH.fij) . The operator }{io'o+d,ai+d,a2+d,a3+d) 
preserve the d + 1- dimensional space V^ao-d-ai-d-a2-d,-a3-d, because — ctj — d E 
{-{ai + d),ai + d+1} and - Y.i=oi~'^i ~ '^)/2 = d. 

Proposition 3.4. (i) We have 
(3.7) 

TT{lo,h,h,h) T — T TT{ao+d,ai+d,a2+d,a:i+d) 

-'-'—ao—d,—ai—d,—a2~d,—a3—d qq— d,— ai— d,— a2— d,— 03— d 

(ii) The characteristic polynomial of the operator H^'-"''-^''-^'^^^ on the space Vag,ai,a2,a3 
coincides with that of the operator 11^°'°^'^'°'^^'^'°''^^''''°'^^'^'^ on the space V-a^j-d-ai-d-a2-d.- 

Proof. Set li = ai + d {i = 0,1,2,3). Then — XlLo^*/^ ~ ^- Theorem 13. 3t we 
have 

jj{-io+d,-li+d,-l2+d,-i3+d) _ _ _ - =2^ - - - - Jj{h,li,h,l3) 

—lo,—li,—l2,—h —lo,—h,—h,—h 

Hence we obtain Eq. (j3.7p . We will prove (ii) in the appendix. □ 

Let fi{x, E) and /2(a;, E) be a basis of solutions to the differential equation i^H^^o,h,h,h) . 
E)f{x) = 0. Since the operator f/'('0''i'^2,i3) jg doubly-periodic, the functions /i(x + 
2iUk, E) and f2{x + 2ujk, E) {k = 1, 3) are also solutions to {H^^oMMM _ E)f{x) = 0, 
and they are written as linear combinations of fi{x, E) and f2{x, E). Hence we have 
monodromy matrices 



(3.8) (A(x + 2u;k, E) f^ix + 2u;,, E)) = (/i(x, E) f^ix, E)) 



(k) (fc) 

'^11 ^12 

(k) (k) 

0-21 '^22 



Now assume that ao, ■ ■ ■ , as, d satisfy the assumption of Theorem 13.31 Set fi{x, E) = 
Lao,ai,a2,a3fi{x,E) {i = 1,2). It follows from Eq. ()3.(ip that fi{x,E) and f2{x,E) are 
eigenfunctions of the operator ("0+^,01+^,02+^,03+^) -^^ith the eigenvalue E. If E is 
not an eigenvalue of H'^^o,h,i2,h) ^j^e space V^o,ai,a2,«35 then {/i(x, E)^ f2ix, E)} is a 
basis of solutions to the differential equation (if {°o+d,«i+d,"2+d,a3+d) _ E)f{x) = 0. It 
is shown in Proposition l3.2l that operator LaQ,ai,a2,a3 is doubly-periodic, and it follows 
from Eq.dnSI) that 

/ (fe) (fe) 

(A(x + 2^,, E) /2(x + 2u;,, E)) = {f\{x, E) E)) 

Hence the monodromy structure of if^''^''^''^'^^^ coincides with the one of f/'(°o+d,ai+d,a2+d,a3 
for the case that E is not an eigenvalue of H'^^o,h,i2,h) ^j^e space VaQ,ai,a2,oi3- If E is 
an eigenvalue of H'^^o,h,h,h) i^j^e space V^o,ai,a2,a3 5 it is also the eigenvalue of the op- 
erator }j{(^0+d,ai+d,a2+d,a3+d) ^^iq spaCe V-ao-d-ai-d-oi2-d-a3-d, which folloWS from 
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Proposition l3.4l fii). Thus the operator Lap, ,02,03 defines an isomonodromic transfor- 
mation from H^^0,h,l2,h) to jjiao+d,a,+d,a2+d,»Wd) _ Similarly L_ao^d,-a^-d,-a2-d,-a3-d 

defines an isomonodromic transformation from ^/'("o+'^.ai+'^.aa+rf.as+d) }j{io,h,h,h) _ 

On the multiplicity of eigenvalues of the operator H^^o,h,h,h) the space VaQ^ai,a2,ai^ 
we have the following proposition: 

Proposition 3.5. (c.f. Let ai be a number such that ai = —U or ai = U + 1 

for each i e {0, 1, 2, 3}. Set d = — Y^i=o o^nd assume d G Z>o and 7^ aj for 
some i,j G {0,1,2,3}. Then zeros of the characteristic polynomial of the operator 
}jih,h,hh) on the space Vao,ai,a2,a3 are distinct for generic periods (2^011,2^73). 

We prove this proposition in the appendix. Now we give a remark on the meaning 
of "generic". In the appendix, it is shown that there exists a periods {2ui,2u3) such 
that zeros of the characteristic polynomial are distinct. Hence the discriminant of 
the characteristic polynomial is not identically zero, and the set of periods such that 
the discriminant of the characteristic polynomial is not equal to zero is open-dense. 

4. Integer case 

We investigate quasi-solvability and generalized Darboux transformation for the 
case k G Z>o {i = 0, 1, 2, 3). Throughout this section, assume U G Z>o for i = 0,1, 2, 3. 

Let be the space spanned by meromorphic doubly periodic functions up to signs, 
namely 

<:i,<:3=±l 

J^ei,e3 = {f{x): meromorphic \ f{x + 2uJi) = eif{x), f{x + 2uj^) = €3/(0;)}, 

If G Z (2 = 0, 1, 2, 3) and — J2^=o '^^/^ ^ ^>0) then Vao,ai,a2,a3 is a subspace of J-'ei,e3 
for suitable J-'ei,e3 (cijCs ^ {il}); because (p(x) — CiY^"^ = pi{x) {i = 1,2,3), where 
pi{x) is the co-p function. Invariant subspaces of with respect to the operator 
HihMMh) is studied in (see also 0011]). 
Let aie{-li,li + l} {i = 0,1,2,3), 

{0}, otherwise, 

and 

{Lao,ai,a2,a3: X^i=0 '^^/^ ^ Z<o, 
— ao,l— "1,1— 02,1— «3 ' 

1, otherwise. 

If Iq + h + h + I3 is even, then the operator ij('o,'i,'2,i3) (^see Eq. p.3|) ) preserves the 
spaces 

Each space is contained in J-'e^,e3 for some ei,e^ G {±1}, and the correspondence 
between the spaces and the signs of (ei, 63) is one-to-one. Let V be the sum of these 
spaces. Then V is written as the direct sum of these spaces, i.e. 

(4.4) V = U-l^-l^-l^-l^ © U-Iq -h,h + l,h + l ® f^-/o,«l+l,-/2,«3 + l ® f^-«0,«l+l,/2 + l,-«3' 



(4.5) g 
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and it is the maximal finite-dimensional invariant subspace in with respect to 
the action of the operator _f/'('o,^i,^2,^3)_ ^^j-^g rearrangement of U such that 

ko > ki > k2 > ks{> 0). If ko + ks > ki + k2 (resp. ko + k^ < ki + ^2), then the 
dimension of the space V is equal to 2/co + 1 (resp. fco + ^1 + ^2 — ^3 + !)• Set 

fco, ko + k3> ki + k2] 

(fco + fci + fc2 - fc3)/2, fco + fcs < fci + fc2- 

Then g G Z>o and dim V = 2g + 1. We set 

(4.6) = (-/o + /l + /2 + /3)/2, = (/o - /l + /2 + /3)/2, 

^2 = (^0 + h-l2 + h)/2, II = {h + I1 + I2- h)/2. 

Note that Iq + + 12 + = Iq + h + h + h ^ 2Z. It follows directly from Proposition 
EH that 

Thus it is shown that the operators which are linked by generalized Darboux trans- 
formations from if ('o''i''2'^3) g^j^g 

(4_7) i/Ca.'f.'i.'i)^ ^{^^2)^ HmUhil)^ H^irhif^)^ 

j^{lo,h,h,h) Jj{h,lo,h,h) jjih,h,l-o,h) j^{h,h,h,lo) 

As is discussed in section |3J these eight operators are isomonodromic. Note that the 
operators ii'('i''0''3''2)^ }j{h,h,io,h) ^ u{h,i2,h,io) obtained by the shift x x + uji 
(i = 1, 2, 3) from the operator H^io,h,i2,h) _ 

Assume that /o, ^1, ^2, ^3 G ^, ^0 > ^1 > ^2 > ^3 > and Iq + li + I2 + I3 is even. 
Set Iq = /g, li = I2, h = li and I3 = max(/Q, — /q — 1). Then /o > ^1 > ^2 > ^3 > 0, and 
the operator H^^o,h,i2,h) jg isomonodromic to iJ('o.'i,'2,«3)_ Note that, if lo + h li+ h 
(resp. jo + ^3 = ^1 + ^2), then we have (/o, ^i, h) 7^ (^~o, h, h) (resp. {Iq, h, k, h) = 
(/o, ^1, h, h))- 

We consider the case lo + h + h + h- odd. If I0 + I1 + I2 + h is odd, then the operator 
]^ik),h,i2,h) prgserves the spaces 

(4.8) ^-«0 -il -i2,i3 + l' ^-io,-il,«2 + l -«3) f^-io,h+l -i2 -i3) ^«0 + l ,-^2 ,-^3 ' 

Each space is contained in J-'ei,t3 for some €1,63 G {±1}, and the correspondence 
between the spaces and the signs of (ei, €3) is one-to-one. Let V be the sum of these 
spaces. Then V is written as 

(4.9) V = f/_i(, _i2,«3 + l © f^-«0 -il,«2 + l -/3 ® + l -/2 -i3 ® Ul^^ + l,~h,-l2,-h^ 

and it is the maximal finite-dimensional invariant subspace in JF with respect to 
the action of the operator ii'('o>'i>'2,'3)_ Lg^ j^. ^^j-^g rearrangement of k such that 
ko> ki> k2> fc3(> 0). If fco > fci ^2 + ^3 + 1 (resp. fco < fci + fc2 + fc3 + 1), then 
the dimension of the space V is equal to 2fco + 1 (resp. fco + fci + ^2 + ^3 + 2). Set 

. / ko, ko>ki + k2 + k3 + 1; 

^ ' ^ \ (fco + fci + fc2 + fc3 + l)/2, fco < fcl + fc2 + fc3 + l. 
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Then g G Z>o and dim V = 2g + 1. We set 

(4.11) /° = (/o + h + l2 + l3 + l)/2, = (^0 + h-l2-h- l)/2, 
^2 = (^0 -h + k-h- l)/2, = (^0 -h-h + h- l)/2. 

It follows from Proposition 13.61 that 

It is shown that the operators which are linked by generalized Darboux transforma- 
tions from if('o:'i.'2,;3) Q^j-e 

(4.12) H^W°2^m^ ffiWlii)^ h(i°2,WD^ H^Wl,i°o)^ 

j^{lo,h,h,h) j^{h,lo,l3,h) j^(l2,h,h,h) j^{h,h,h,lo) 

and these eight operators are isomonodromic. 

Assume that Iq, li, I2, ^3 G Z, /q > ^1 > ^2 > ^3 > and Zq + ^1 + ^2 + ^3 is 
odd. Set !o = 1°, h = 1°, I2 = l°2 and [3 = max(/°, -/° - 1). Then [q > h > h > 
h > 0, and the operator ^^"(^0,^1,^2,^3) jg isomonodromic to if('0''i''2,'3)_ Note that, if 
^0 7^ ^1 + ^2 + ^3 + 1 (resp. jo = h + l2 + h + ^), then we have (Zq, h, k, h) 7^ (^~o, h, h, h) 
(resp. {lo,h,l2,h) = {lo,h,k,h))- 

Now we reproduce results in this section for the case of Lame equation and asso- 
ciated Lame equation. If three (resp. two) of k {i = 0, 1, 2, 3) are zero, then Eq. ()1.4|) 
is called the Lame equation (resp. the associated Lame equation). For simplicity, we 
consider the case I2 = h = 0, /q, /i G Z and Iq>Ii> 0. 

If Iq + li is even and /q > ^i, then the operator if('o,'i,o,o) jg jgomonodromic to 
^((«o+ii)/2,(Zo+«i)/2,(«o-«i)/2,(/o-/i)/2-i) ^Yie transformation L_,o,,^+i,i,o. Especially, if 
Iq is even and /i = (the case of Lame equation), then ffC^o, 0,0,0) isomonodromic 
to if('o/2,«o/2,«o/2,«o/2-i)_ j^Qi-g ^^g^i-^ jf jg gygj^ ^ the operator 

ffiioMm jg gelf-dual. 

If Iq + li is odd and Iq > li — 1, then the operator 7f ('o,'i,o,o) jg jgomonodromic 
to ff{(io+h+i)/2,ao+/i-i)/2,(/o-/i-i)/2,{io-h-i)/2) ^j^g transformation L_io,i,+i,i,i. Espe- 
cially, if Zo is odd and Zi = (the case of Lame equation), then ^^('0,0,0,0) jg jgomon- 
odromic to i7(('o+i)/2,(Zo-i)/2,(io-i)/2,(/o-i)/2)_ Note that, if /q + h is odd and lo = h- 1, 
then the operator Co, '1,0,0) jg gelf-dual. 

Thus we confirmed the conjecture of Khare and Sukhatme |SI on isospectral partner 
of Lame equation and associated Lame equation. 

5. Half-integer case 

We investigate quasi-solvability and generalized Darboux transformation for the 
case /, + 1/2 G Z (i = 0,1,2,3). Write k = rii — 1/2 and assume G Z>o for 
i = 0,1,2,3. 

Candidates for the space related to quasi-solvability written as Eq. ()3.1|) are de- 
scribed as l^g+i fj^+i fj^+i.ns+ij where hi G {±ni} {i = 0,1,2,3). For the existence 
of non-zero space V^^+i j^^+i j^^+i.ns+i) the condition X]i=o^« ^ necessary. In 
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other words, if ^j=Q rii is odd, then the space related to quasi-solvabihty written as 
Eq. ljH.lj) does not exist. 

On the rest of this section, we assume that XlLo'^* even. Then the operator 

^("■0-2 '"1-2 '"-2- 2 '"-3- 2) preserves the spaces 



-no+5,-ni + i,-n2+5,-n3+5 ? ^-no + 5,-"-i + 5>"2 + 5,n3+^ ' 
-".0+|,"i + |,-i.2 + f ,"3+1 ' ^-no+|,rai + |,n2 + |,-n3+| ) 

TT 

-no+5,-ni + ^,-n2+5,n3+^ ' ^ -no+^,-ni + ^,n2 + ^,-n3+-^ ' 

TT 

-no+5,ni + i,-n2 + |,-n3+i ' no+|,-ni + i,-n2 + |,-n3+i ) 



where Uao,ai,a2,a3 is defined in Eq. (l4.1|) . Unhke the case lo,li,l2j3 £ ^5 these 
eight spaces are not disjoint. In fact, these eight spaces are subspaces of the space 



^0 = (no + ni + n2 + ^3)72, 
= {no -ni + n2- n3)/2, 
= {-Uq + ni + 77-2 + ^3)72 

n\ — {riQ + ni- n2 + ns)/2, 



nl ' = {no + ni-n2- ^3)72, 

^3^^ = ('^O - ni - 722 + n3)/2, 

nf^ = {no ~ni + n2 + n^)/2, 
"^3 = {nQ + ni+n2- n^)/^. 



It follows from Proposition 13.61 that 



(1) 



^3 



"0+0 



-n2^ 



(1) 



1 „(!)- 
2' 2' 3 



1_ 

" 2 ' 

2 ' 2 



-"3 + 



"0 + ^, 
1 



-n2- 



^>"2- 



'"3 + ^ 



;,"3- 



■'^0+2 >"l+2 '-"2 + T 



-"3 + ^ 



(1) 



-"0+2 
2 '"3 1 



-"2 + ^, -"3 + 



liy("0-5,«i-|,n2-|,n3-|) 



^0 



^-no+i, 

(1) 



-ni + |,n2 + |,-n3 + ^ 



no+^,-ni + 5,"2 + |,-n3+| 
ni-5,n2-5,n3-5) 



^0 



-"0 + 7.-"l + i,-n2+^,»i3+^ 



-"■0+9 ,-ni + T 



1 



r."3- 
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where L 



is defined in Eg. ()4.2|1 . We also obtain that 



(2) 



,(2) 1 „(2) 1 „(2) In ~ 



-"0- 



2'"'0 2 '"-3 2 '"-2 2''/; 



-"■0 + 5,-"l + 5,-'"2 + 5,-ri3+5 
lij("0-i"l-|>"2-i«3-|) 



-Jii + i n2 + 4,n3+ 



■rao+^,-ni + i,n2 + |,n3+i 
5,n2-5,n3-|) 



/ (2) 
_ff("2 - 

= L_ 
= L 



,(2)_1 „(2)_ 



-no+|,rii + |, 
1 „(2)_1 „(2)_ 



,(2) 1.^ 



-112- 



r,"2- 



m"3- 



I) 



,(2) In; 



'0 



-n3- 



-713- 



iif("0-i"i-5'"2-5,n3-5)_ 



Thus it is shown that the operators which are hnked by generahzed Darboux trans- 
formations from i7("o-2'"'i~2'"2-2>"3-2) are 



^(„W„i,„(i)_i,„(i)_i,„(i)_i 

rrfn*^)- 1 n<^'-i n<^'-i n<^'-i 
_f^i"2 2 '"-3 2'"l 2'"'0 2 

_f^{n2-5,n3-5,no-i,ni-i) ^ 



2 '"O 2 '"3 



2 '"2 



(2)_1 „(2)_ 



,(2)_1 „(2)_ 



(2) 1 (2) 



rai-|,no-|,n3-|,n2-|) 
™3-5,n2-5,ni-5,no-5) 



{2)_i (2) 



Now we investigate the case Uq, rii, n2, ns G Z>o, uq + rii + n2 + e 2Z>o and 
no > Ui > n2 > n^. Recall that li = Ui — 1/2 {i = 0, 1,2,3). We divide into three 
cases, the case uq > rii + n2 + n^, the case rii + n2 — < rio < rii + n2 + and the 
case Uq < Til + n2 — ""-s- if ao + c^i + «2 + «3 = 2, then we set Vao^ai,a2,a3 = {0}. 
If no > ni + n2 + n^, then the operator if ("o^l'"i^l'"2-2."3-|) preserves the spaces 



Vi 

2 



-no i-ni i-?i2,|+n.3 



-"2, 



'"3 



h-no,i-ni,^+n2,i+n3 




-?i2,i+n3 



9 -no, 



2 no,|-|-ni,|+n2,|-?i3 



Here V. 



ao,ai,a2,a3 



^"^ctQ CXI CX2 CK3 means that ^''^q (^\.c^2-c^3 ^ subspace of ^^^q^o qi ck2 ^^3* On 



this case, the operator if ("o 2 '"-i 2 '"2 2 2 ) is isomonodromic to ii'^^o 



,{2) 



and H^^3 2 '"2 2'"'i 2' 



'0 2 



by the transformations L_ 



no 



+ ^,ni-|-i n2-|-4,n3-|-i 



and 



-^-no-|-|,ni-|-|,n2-|-|, 

If Ui + n2 — ^ no < Ui + n2 + n^, then the operator ii("-o-2'"i~2'"2-2."3-2) 
preserves the spaces 
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^i-no,^-ni,i-n2,i-n3 



Vl 



9 -no 



,^-ni,^+n2,i+n3 



2 no,5+ni,5-n-2 



"3 



-"0,5 +"1,5 -"2, 5+^3 



V: 



i-no,4+ni,^+n2,^-n3 ' 



Vl 



+no,:T-ni, 



and it IS isomonodromic to -H 2'"i 2 '"2 2 '"3 2''andiJ'> 3 2 '"2 2'"i 2'"o 2 J 



by the transformations __„,+ i 

If no < rii + ^2 — ris, then the operator if 

Vl 



2 "u ' 2 '"-^ ' 2 

no-^,ni-5,n2-5,n3-|) 



1^1 



-no, 



-711, 



-"•2, 



-713 



_ -no, I -ni, I -712, 1 +"3 
^i-no,|-ni,5+7i2,5-n3 
^5-no,|+ni,|-n2,|-n3 
^^+no,^-ni,^-n2,5-n3 




preserves the spaces 



i-7io,i-ni,i+n2,i-|-n3 



- 710 , i ~(-ni , i - 712 , i -|-n3 



§+7io,§-ni,|-»i.2,|-l-n3 ' 



and it is isomonodromic to if("o 
by the transformations L^^ 



2'"l 2 '"2 2' " 



(1) In 
3 -2l 



-"3- 



) and if("3 
1 and L 



(2) 1 m 1 ^(2) 1 ^(2) 1 



-"6 



-n2-f ^,"3 + ^ ■ 



6. Finite-gap potential 

In this section, we construct an odd-order differential operator A (see below) which 
commutes with the operator if ('o''i''2,'3) composing generalized Darboux transfor- 
mations for the case lo,h,l2,h G ^>o- We also show that the commuting operator 
coincides with the one constructed in 18 . We use notations in section IH 



-^-«|,/^-(-l,«f + l,-«,'j-^-«l,/o + l,-/3,/2-Hl-^-«0,-if,/i-hl,/|-(-l-^-«O,-/l,-/2,-i3- 



-^«g+i,-/g,-«S,-«j-^-h,-«o,i3+i -^2-^-^(';,/^+l,-^g,-^g-^^o-hl,-/l,-^2,-^3• 



Proposition 6.1. If Iq + h + h + h is even, then we set 

(6.1) A 

If lo + h + h + h is odd, then we set 

(6.2) A 
The operator A commutes with _f/'('0''i''2,^3)^ ^ 

Proof. If Iq + h + I2 + h is even, then Eq. ()6.3p is shown by applying the following 
relations: 

(6.4) ii('§''f''i''i)Z.,„,_,„„,,,_,3 = L_,„_,,,_,,,_,3ii('«''-'^''«), 

ij ^ ' ' ° i--/i,io+l,-/3,i2+l - -^-iiAi+l -i3,/2+l-n' S 

^^'°''^''^''^^^-;|,;^+i,zf+i,-;g =^-z|,/i+i,/f+i,-ze^^'^''^''^''«^- 
If Iq + li + I2 + I3 is odd, then Eq. (j6.3|) is shown by 

(6.5) ii('0°,'?,'2,'3)L,,^,,„,,,_,,,_,3 = L,, + l,_,„_,„_,3^('-'-'-'3)^ 



jj{h,lo,h,h)l ^ I 

'oi'lT^J-, i2' '3 'o' 1^ ' 



^{l0,h,h,h)T 

to 



Zj{lZ,l'l,l°2,ll) 



■h,—lo,l-i+'^,—h 



■'l, — '0,'3 + l, — ^2 



lj{h, lo,h,h) 
,-,±1 , 



2+l'~'3'~'0'~'l 
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□ 

Note that the order of the differential operator A is equal to 2(7 + 1, where g is 
defined in Eq. (l4.5j) or Eq. ()4.1()j] . Recall that the space V is written as Eq. ()4.4|l or 
Eq. ()4.9j) , dim V = 2g + 1 and it is the maximal finite-dimensional invariant subspace 
in with respect to the action of the operator ff('0''i''2>^3)_ 

Proposition 6.2. The operator A is the monic differential operator of minimum 
order which vanishes all elements in V . 

Proof. We show this proposition for the case that I0 + I1+I2 + h is even. For the case 
that Iq + h + I2 + h is odd, it is shown similarly. 

Since V is written as Eq. ()4.4|) and degA = dimV^, it is enough to show that 

The operator A is written as Eq. ()6.H) . and the operator L^i^^^^i^^^^i^^^i^ vanishes any ele- 
ment in the space U-if^-i-^-i^-i^. Hence we have A(f){x) = for G U-i^-i-^-i^-i^. 
Let {fi{x), . . . , fn{x)) be a basis of the space -/^^jj+i^^g+i. Since the operator 
}j{k)AMM preserve the space U^i^^_i^^i.2+i,h+^-: function H^^°'^^''-'^'^^'> fj{x) is written 
as Yl^=i (^i,jfi{^) for some constants a^j-. It follows from Eq. ()6.4|) that 

n 

= (^i,jL-la-l^-l2-lJj{x). 
i=l 

Set U,e = -^-io-h -«2 -«3f^-fe -;i,;2+i,«3+i- Then it follows from Eq.()6.6|) 

that the space U_i^_ie^ie^-i^ ie^i is invariant under the action of if^'o.^!''!''!). Let (£1,63) 
(ei,e3 G {±1}) be the numbers such that f/--io,-h,«2+i,i3+i C J^ei,e3- Then we have 
ie^_^_i^ie_^_i c J^ei,e:i and U-ie_ie^_^_i^ie_^_i c J^ei,e:i- As is showu in P Theorem 
3.1], the space f/.^e _;e /e_,_]^_;e_,_i is the maximum subspace of J-'ei,e3 which is invariant 
under the action of iir('o''i''2''3). Thus we have 

f^-ig -if,Z| + l,i§ + l = -^-io-h -«2 -«3f^-i0,-«l,i2 + l,«3 + l + + 

Similarly it is shown that 

-^-ig ,-«J,ii + l,«| + l-^-i0,-«l -«2 ,-«3f^-«0,«l + l,-«2,«3 + l <^ + 1 -«3,«2 + l5 

+ 1 -i3,«2 + l-^-ifj,-«f,«i + l,«i + l-^-i0,-h -«2,-«3^-«0,h+l,i2+l -«3 + + • 

Since the operator ^e+i ;e+i (resp. L^i^^i^+i^^i^^i^+i, L.^e je+i ^e+i vanishes 

any element in the space f/.^g ^e+i (resp. U-i, ^1^+1,^13,12+1, f/_/|,/|+i,;j+i we 
have A(j){x) = for (j){x) G U-i^ -1^,1^+1,1^+1 (resp. G U-ig,i^+i,-i^,i^+i, G 

f^-/0,il + l,«2 + l -«3)- 1^ 



It follows from Proposition 16. 21 that the kernel of the operator A coincides with the 
space V. We denote the monic characteristic polynomial of the operator ff('0''i''2.'3) 
on the space V by P{E). For simplicity, we set u{x) = Yl^=o + ^)pi^ + ^i) 
H = HiioMMM = -d^/dx^ + u{x). 

Proposition 6.3. Set do{x) = 1 and dg+i{x) = 0. The operator A is written as 

9 

-j 



(6.7) A 
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for even doubly-periodic functions aj{x) {j = 1, . . . , g) and constants Cj {j = 0, . . . , g), 
and dj{x) {j = 0, . . . , g) satisfies 

(6.8) aj'{x) — 4:u{x)aj{x) + 4:aj^^{x) — 2u{x)dj{x) = 0. 

Proof. Since A is a monic differential operator of order 2(yf + 1, it is written as 



A = {-iy 



E 

.3=0 



wliere ao{x) = 1. We have 



j=0 



JJ9-J 



J2 [-aA^Wix) + mx)^ + i~a';ix) + m^))^ + b';ix) ) H^-^ 
j=o \ X X 

9 



d 



J2[2~a'^ix)i-H + uix)) + i~a';ix) + 2^ix))- , , 

j=0 ^ 



d,ix)u'ix) + b'Ux) I if 

' dx 

J2 ((a'^ix) + 2&;(a;))^ - 2fi;+i(x) + 2d'^{x)uix) + dj{x)u'{x) + b'^ix)^ H^-^ 



j=0 

= 0. 

Hence we have 

aj(x) + 26;. (x) = 0, -25;.+i(x) + 2d'^{x)uix) + djix)u'ix) + b"^{x) = 0. 
Therefore 

bj{x) = —d'j{x)/2 + Cj, d'j'^x) — 4:u{x)d'j{x) + 4:aj_^_i{x) — 2u{x)dj{x) = 
for some constants cj {j = 0, . . . , g), and we obtain the proposition. □ 
Set 



i=0 



(6.9) E{x,E) = Y,a9~iix)E\ 

It follows from Eq. ()6.8p that 
(6.10) 



Set 



— - 4 M X -E) — - 2m' X H X, E = 0. 
ax'^ ax 



1^. rf2H(x,E) 1 / dE{x,E) 



(6.11) Q(^) = S(x, E)' (^ - u{x)) + -S(x, E) 



dx'^ 4 \ (ix 



It is shown by differentiating the right hand side of Eq. ()6.1ip and applying Eq. ()6.10|) 
that Q{E) is independent of x. Q{E) is a monic polynomial in E of degree 2g + 
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1, which follows from the expression for E') given by Eg. Similarly to 

Proposition 3.2 in [Hj, we can show 

Let us recall the function defined in It satisfies Eg. (jfi.lOj) and has an 

expression 

3 h-l 

(6.12) E) = co{E) bf^EMx + u:.f-\ 

where the coefficients cq{E) and lf^\E) are polynomials in E, they do not have com- 
mon divisors and the polynomial Cq{E) is monic. It is shown that the dimension of the 
functions which are doubly-periodic and satisfy Eg. ()6.im) is one (see |jlO,j Proposition 
3.9]). Since H(x, E) is a polynomial with respect to the variable E and coefficients of 
H(x, i?) are coprime, we have 

(6.13) S(x, E) = S(x, E) [y^ diE"' 

\i=0 

for non-negative integer k and constants di {i = 0, . . . , k) such that do = 1. Write 

g-k 

(6.14) E{x,E)=J2a9-k-^{x)E\ 
Then we have aQ{x) = 1 and 

k 

i=o 

The functions aj(x) {i = 0, . . . , g — k) also satisfy Eg. ()6.8|) . because S(x, E) satisfies 
Ea. ljFmH) . Set 

(6.15) ^ = g{„^(.)iL_l(i^„^(,))}H»-- 

It follows from Eg.dHSI) for ai{x) that [A, H] = 0. 
Proposition 6.4. 

9 / k 

(6.16) {-lyA -J2^9-iH' = A [J2dk~jH^ 

1=0 \j=0 
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Proof. 

(-l)M - ± c,^,H' = ± {a.(x) A _ 1 {±aA.)) ] 

1=0 i=0 ^ ^ ^ ^ 



, dx 2 \ dx 

i=0 {. \j=0 / \ j=Q 

g k 



H 



9-« 



f d 1 / d \1 

i=0 j=0 ^ ^ ■' 



vj=0 



□ 



Set 



2,^ 1_ c/22(x,E) 1 fdE{x,E) 



(6.17) g(E) = S(x,i5;)^(E-w(x)) + -H(x,E) 

Then the right hand side of Eq. ()6.17p is independent of x, and Q{E) is a monic 
polynomial in E of degree 2{g — A;) + 1. It is shown in (HI Proposition 3.2] that 
= Q{H). By Eq.dniSl) and definitions of Q{E) and Q{E), we have 

(6.18) g(E) = g(^) l^^t/.E'^-^j , degg(E)<degg(E) = 2(7 + 1. 

On zeros of Q{E) and P{E), the following proposition is shown in [HI Theorem 3.8] 
(see also jHl Proposition 2.4]): 

Proposition 6.5. (c.f. ^ Theorem 3.8]^ The set of zeros of Q{E) coincides with 
the set of zeros of P{E). 

Proposition 6.6. (c.f. [TT]^ Roots of the equation P{E) = are distinct for generic 
periods {2lji, 2uj3). 

We prove Proposition 16.61 in the appendix. 

Proposition 6.7. For the periods {2uji, 2003) that roots of the equation P{E) = are 
distinct, we have P{E) = Q{E) = Q{E) and E{x,E) = E{x,E). 

Proof. By assumption, the equation P{E) = has 2g + 1 roots which are distinct. 
It follows from Proposition 16.51 that the number of roots of the equation Q{E) = 
is equal to or more than 2g + 1. Then we have deg g(i?) > 2g + 1. Combining with 
Ea. dFTTKI) and degP(£;) =2g + l, we have Q{E) = Q{E) = P{E), and the value k in 
Eq. (j6.13| ) is zero. Hence E{x,E) = E{x,E). □ 

Proposition 6.8. All constants Cj {j = 1, . . . , g) are zero. Namely we have 

(6.19) i = (-l)'i:{a,W^-i(^a,M)}i/" 
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Proof. First, we assume that roots of the equation P{E) = are distinct. Let 
{i?j}j=i^..._2g+i be the roots of the equation P{E) = 0. Then Ei {i = 1, ... ,2g + 1) is 
an eigenvalues of the operator H on the space V. Let G be the eigenfunc- 

tion of the operator H with the eigenvalue Ei. It follows from Proposition lfj.4l that 
ii-^y^-Elo^a-iH'? = A' = Q{H) = P{H). Hence (i-2(-l)^' Ef=o + 
iYlf=o ^9-1^^)"^ = P{H). We apply this operator to (pi^x). Since the operator A an- 
nihilates the space V, we have (XlLo ~ P{Ei)(f)i{x) = 0. Hence 
Ylf=o ^a-ii^^y = for 2 = 1, . . . , (7. Since the polynomial of degree less that 5^ + 1 
cannot have 2g + 1 zeros, we have Q = {i = 0, . . . , g) for the case that roots of the 
equation P{E) = are distinct. 

Now consider the case that roots of the equation P{E) = are not distinct. Since 
the operator A is defined by Eq. ljfi.lj) or Eq. ()fi.2|l . the coefficients q {i = 0, . . . , g) are 
continuous as a function of periods {2uJi, 2ujs). Since q = (2 = 0, . . . , (7) for a dense 
set of periods, we have q = (i = 0, ...,(?) for all periods. □ 

For the case that roots of the equation P{E) = are distinct, it is shown that 
P{E) = Q{E), the value k in Eq. (j6T3|l is zero, E{x, E) = E{x, E) and A = (-l)^A. 

We now consider the case that roots of the equation P{E) = are not distinct. It 
is already shown that P{E) = Q{E) for a dense set of periods, deg P(ii^) = degQ{E) 
for all periods, and coefficients of P{E) and Q{E) are continuous with respect to 
the periods. Hence we have P{E) = Q{E) for all periods. Assume that the value k 
in Eq. ()6.13p is positive. Combining with Eq. ()6.18|) and Proposition 16.51 all zeros of 

(X]i=o ^i^^Y s^re zeros of Q{E). Let Eq be a zero of J2i=o d'iE\ i.e. ^j=o dii^oY = 0- 
By Eq. ()6.16|) and Proposition 16.81 it is shown that, if /(x) satisfies Hf{x) = Eof{x), 
then Af{x) = 0. Hence all solutions to {H — EQ)f{x) = are contained in the space 
V. But it contradicts that all solutions to {H — EQ)f{x) = cannot be doubly- 
periodic up to signs, which is essentially shown in [6^ Theorem 3.8] (see also the 
proof of [TUl Proposition 3.9]). Hence we have k = 0. Therefore A = {—1)^A, 
deg P{E) = degQ{E) and P{E) = Q{E) for all periods. It follows from Proposition 
16.21 that the operator {—iyA{= A) is characterized by the monic operator of order 
2g + 1 which annihilates all elements in the space V. Summarizing, we obtain the 
following theorem: 

Theorem 6.9. Let A be the operator defined by composing generalized Darboux trans- 
formation (see Eq. \6.1\) or Eq. \6. jj)) ) and A be the operator defined from the even 
doubly-periodic function E{x,E) (see Eqs. l[6.1^ W.15}) ). Let P{E) be the characteris- 
tic polynomial of the operator H on the space V , and Q{E) be the polynomial defined 
from E{x,E) (see Eq. \6.11\l ). 



(i) We have A = (-l)M and P{E) = Q{E). 

(a) The operator {—1)^A is also characterized by the monic operator of order 2g + 1 
which annihilates all elements in the space V. 

Note that we proved Conjecture 1 in jSj by (i) and Conjecture 2 in |5] by (ii). Since 
P{E) = Q{E), the genus of the spectral curve z/^ = —Q{E) is g, where g is defined 
in Ea. dlB (/q + + /2 + ^3: even) or Ea. (ICT]|l (/q + h + h + h: odd), and it agrees 
with the results in [21 ^] . 

We consider isomonodromic property again. Let , Z° = 0, 1, 2, 3) be the numbers 
defined in Eqs. ()4.61 14. lip . We have shown in section HI that, if Iq + li + I2 + I3 is 
even (resp. odd), then the operator if ('0''i''2''3) jg linked to _f/'('0''i''2>'3) generalized 
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Darboux transformation, where //('o-'i'^a.^a) jg g^^y operator listed in Eq. ()4.7p (resp. 
Eq. ()4.12j) ). We now show that functions related to monodromy for the operator 
jjik),h,h,h) coincides to ones for the operator if ('o''i''2.'3)_ Pqj^ ^]-^jg purpose, we recall 
functions defined in IH1 E] ■ 

Let S(^0'^i''2''3)(a;, E) be the function defined in Eg. fjO^ . Q(^o,iuh,i3)(^E) be the poly- 
nomial defined in Eq. ()6.17|) and p('-0''-^'^^M(^E) be the polynomial P{E) defined in this 
section. Set 



^{io,h,i2M[x,E) 

Then the function A(^oM2M(^x, E) is a solution to the differential equation 

^ffiioAMM _ E)f{x) = 0. 

Set 

a) = exp(C(a)x), = 0, 1, 2, 3), 

where a{x) is the Weierstrass sigma function. The function A^'-°''-'^''-'^'''^\x, E) admits 
an expression in terms of Hermite-Krichever Ansatz. Namely, we have the following 
proposition: 

Proposition 6.10. (c.f. Theorem 2.3]^ The function A^^'^''^^''''^'^^\x, E) is expressed 
as 

(6.20) A('''''i''2''3)(x,E) = exp(fi:x) (j-)' 

for some values b^j^ (z = 0, . . . , 3, j = 0, . . . , Zj — 1) except for finitely-many eigenvalues 
E. The values a and k are expressed as 



p(a) = R^^''''''''''\e), p'{a) = R^^°'''''''''\e)^-Q('oM2M(e), 

K = R^^'"^"^''^''\e)^-Q(^oM2M(e), 

where 4'°'''''''''^(^), 4'°'^'''''''^(^), 4'"''''^''''^^) are rational functions m E. 
It follows from Eq. fjFT^ that 

(6.21) A('oA,i2,«3)(3, ^ 2uJk, E) = exp{-2rika + 2uJkC{a) + 2KuJk)A^-^'''^''^''^'\x, E) 

for = 1, 2, 3. A monodromy formula in terms of hyperelliptic integral was obtained 
in [H|. 

Proposition 6.11. (c.f. |Hl Theorem 3.7]^ Let k G {1,2,3}, qk E {0,1} and Eq be 

the eigenvalue such that hiio,h,i2,h)(^x + 2uk, Eq) = [-l)'ikA{io,ii,i2,h)(^x, Eq). Then 

(6.22) A'^ioMhM^^ + 2u}k, E) = (-1)^'=A('0'^^''2''«)(x, E)- 

exp -- / / = dE 

\ ^JEo ^_Qilo,li,h,l3)(E) 

where a^^O'^^'^^'^^\E) (resp. c'^^o,h,i2,h)i^E)) is a polynomial defined in [Oj. 
The following proposition states the coincidence of functions. 
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Proposition 6.12. Let ij('o:'i.'2,;3) g^^y operator listed in Eq. d^. ?| j (Iq + I1 + I2 + h: 

even) or Eq. (Iq + h + I2 + h- odd). Then we have 

p(joM2M^^-^ ^ p{hhhk)^p^^^ ( • _ 2, 3), 

Proof. We show the proposition for the case that /o+^i+^2+^3 is even and (/q, /i, ^2, ^3) = 
(Z3, ^1; ^o)- -^o^ other cases, it is shown similarly. 
The space V for the operator /7('0''i''2.'3) jg -^yntten as 

and the corresponding space for the operator if(^3''2''i''o) is written as 

It is seen that the space U^i,^_i^_i^_i.^ (resp. f/.^^ .Zj^^^+i.^s+i' f^-«o,/i+i -/2,/3+i' f^--/o,'i+i,/2+i, 

is linked to the space f/_;|_i|_i|_ig (resp. f/_;e _;e ;e + i ;e + i, f/.^e ,e + i _;e + [/(e + i _;e_^,e ,g_^i 

by the generalized Darboux transformation L_i^^_i^_i^_i^ (resp. L^^^ i2+i,/3+i5 

^-/o,«i+i -i2,/3+i> ^-«o,/i+i,«2+i -is) and the shift x x + (resp. x ^ x + uji, 
X — > X + u;2, X — > x). It follows from Proposition 13.41 that the characteristic polyno- 
mial of H^^oMMM on the space U-i^-i^-i^-i^ (resp. U-i^-i^^i^+i^i^+i, f/_i„,;,+i _i2,/3+i, 
f/_io,i,+i,/2+i ^ig) is equal to that of if(^3.'i.^f.^o) on the space U-^-q-ie-ie (resp. 

f/_ie ;e + i + f/.^e ;e + i „je + i , f/i| + 1 _ ,e _ ;e ^g + i ) . SlUCC the polyUOmlal P {E) IS Writ- 

ten as the product of characteristic polynomials of invariant subspaces, we have 
p{io,h,i2,i3)(^E) = P^Whn:)[E). It follows from Theorem EH that Q('o''i''2''3)(^) = 

Let a and k be the corresponding values in Eq. (|6.21| ) for the parameters (Zg, Zg, Z^, Zq). 
Namely, they satisfy 

KWhil){^x + 2uJk,E) = exp(-2r7fea + 2cUfeC(«) + 2/icUfc)A('3.'2.'?.'o)(a;,^) 

for k = 1,2, 3. The monodromy matrix is preserved by the generalized Darboux 
transformation for almost eigenvalues E. Since the values exp(ib(— 277^0 + 2ujk({a) + 
2KLjf:)) (resp. exp(ib(— 2?7fca + 2LJk({c() + 2kiJk))) are eigenvalues of the monodromy 
matrix, we have 

- 2r]ia + 2cJi(C(a) + k) = ±{-2T]ia + 2cJi(C(a) + k)) + 2nx-n\f-\, 

- 2ri3a + 2cj3(C(a) + k) = ±i-2ri3a + 2cj3(C(5) + k)) + 2^371 v^, 

for some integers ni and and almost eigenvalues E. By the Legendre's relation 
rjiuj^ — ?73u;i = 7rv^^/2 and the relation ({x + 2uk) = C{x) + 2rjk {k = 1,3), it follows 
that 

a = ±a — (2niu;3 — 2n3u;i), k = ±k. 

It follows from the asymptotic of k (see jHl Proposition 3.2]) that the sign ± is plus. 
Hence p{a) = p{a), p'{a) = p'{a) and /t = k for almost E. Since /^('0''i''2''3)j^^-j 
i?f^''^'''''o^(E) {i = 1,2,3) are rational functions, we have /?f°'''''''''^(E) = i^J's.'i.'f.'o)^^) 
for z = 1, 2, 3. 
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Let Eq be the corresponding values in Eg. ()fi.22|) for the parameters l^y ^ly ^o) ■ 
By applying a similar discussion for Eq. ()(i.22|l . we obtain that the integrals 



a 



dE- / — — dF, 



— , ^dE — / — , ^dE, 

Eo J _Q{h,hMz){^E) "^^0 J ~Q^Whio){E) 



are constants. By differentiating them in the variable E and using the relation 

Q{io,h,i2,h)(^E) = g('3''i.'?.'o)(E), wehavea('0''i''2,«3)(;^) = a^Wl^^D^E) and c^^°'^^'^^^^^\E) 
c('i.'i''!''8)(^). □ 

It follows from Proposition 16. 12l that the hyperelliptic-to-elliptic integral formulae 
obtained in fQ^ for the parameters {loJij2,h) coincide with that for the parameters 
(^0, ^1, h, h)- 

We calculate the commuting operator A explicitly for the case lo = h = h = h = 
g E Z>i. The commuting operator A of H^^'^'^'^^ is written as 

r 1 d\ - 

(see Eq.dH), and = - ei)-a/\p{x) - e2)-^/'(p(x) - 63))-^/^ = 

(p'(x)/2)^^. Hence we have 

(6.23) ^ = (-l)VW-o(-J-^iL)^"'o 

On the other hand, we have X]i=o ^'(•^ + ~ 4p(2x). By the change 2x — ^ x, we 
recover the Lame operator 4(— (i^/(ia;^ + (7(5f + l)p(x)). Hence the commuting operator 
A of i7(f '0,0,0) jg written as 

(6.24) A = {-iyp'{x/2y+' o o p'{x/2y. 



p'[x/2) dx J 
Therefore we obtain the following proposition: 

Proposition 6.13. The commuting operator A for the Lame operator iJ(5''0'0,o) (z 
Z>i) is written as Eq. \6.24 )- 



Appendix A. Elliptic functions 

This appendix presents the definitions of and the formulas for the elliptic functions. 
The Weierstrass p-function is defined by 

(A.l) p{x)- — + ^ \{x- 2ma;i - 2^^3)2 ~ {2muJi + 2^^3)2 ) ' 

(m,n)eZxZ\{(0,0)} 6/ \ i 6/ / 



HEUN EQUATION V 21 

Setting ijj2 = —uJi — Us and = p{ujk) {k = 1, 2, 3) yields the relations 
(A.2) ei + 62 + 63 = 0, 

2Vp(a;)-ei p(x)-62 p(x)-63/' 
+ =ei + ^^ ^4^^ — z = 1,2,3, 

where i', i" G {1, 2, 3} with < i" , i 7^ i' and z 7^ i". 

Appendix B. Proofs of Propositions 13.41 (ii). 13.51 and IH.HI 

Let ai be a number such that ai = —U or = /j + 1 for each i G {0, 1, 2, 3}. Set 
a = (ao, «!, 02, as) and 

<+i,, = -4(r + 7f)(r + 72^), 

0^r-\,T = -'irir + a2- l/2)(62 - 63)(62 - 6i), 

a^_^ = -4r((62 - 63) (r + 0:2 + ai) + (62 - 61) (r + 02 + 0:3)) 

- 4627^72 + ei(a2 + as)^ + e2(«i + a^f + 63(01 + a2)^ 

71 = («o + ai + 6^2 + a3)/2, 72 = ("C^o + ai + ^2 + "3 + l)/2. 
Then the action of the operator if ('o.'i:'2,«3) jg written as 
(B.l) H(io,hh,h)y^ = <+iX+i + <X + <-i,r<-i' 

(see j7]). Set c? = — X^Lo'^*/^ ^^"^ assume d G Z>o. Then the space V^o, 0-1,0-2,03 
is spanned by f g , f f , . . . , f J, and it follows from Eg. (jB.lj) that /7('0''i''2:'3) preserves 
the space V^o,ai,a2,a3- The operator }ji"'o+d,ai+d,a2+d,a3+d) g^jgQ preserves the d + 1- 
dimensional space V_ao-d~ai~d,^a2-d-a3-d- 

Proposition B.l. (Proposition \3.4\ (ii)) Set d = —Yl^=o'^i/'^ ^^'^ assume d G 
Z>o. Then the characteristic polynomial of the operator if('0''i''2''3) ffi^ space 
K«o,oi,a2,a3 coincides with that of the operator }j{°'o+d,ai+d,a2+d,a3+d) space 

^—ao—d,—ai—d,—a2—d,—ag—d ■ 

Proof. Set —a — d = {—aQ — d, —ai — d, —a2 — d, —a^ — d). The action of the operator 
^(ao+d,ai+d,a2+d,a3+d) ^-^e space V"_oo-d,-ai-d,-02~d,-03-d is Written as 

TT{ao+d,ai+d,a2+d,a3+d) -a-d _ -a-d -a-d , -a-d -a-d , -a-d -a-d 
11 Uj. — '^r+l r r r ' ^r— 1 r 

Set uf = {{—lYd\/{r\{d — r)\))v^_,^,. By a direct calculation, Eq. ljB.lll is rewritten as 

Hence the martix representation of if ("o+d,ai+d,a2+d,a3+d) the basis {v(^''~'^, • • • , v^''^'^) 
is the transposed one of if Co.'i-'a.^s) ^^j^g basis (mq, . . . and we obtain that the 
characteristic polynomial of the operator ii'('o.'i.'2,«3) ^j-^g space V^o,cn,a2,03 coincides 
with that of the operator ii("o+'i.ai+'^>«2+d,a3+a!) ^^le space ■l/_oo-d,-oi-d,-02-d,-a3-d- 

□ 
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Proposition B.2. ( Proposition \8. ,5j) Set d = — Yl^=o^i/'^ '^'^^ assume d G Z>o and 
ai 7^ aj for some i,j G {0, 1, 2, 3}. Then zeros of the characteristic polynomial of the 
operator H^'""''''^'^^'''^^ on the space V"ao,ai,o2,o3 distinct for generic periods {2uJi, 2UJ3). 

Proof. The assumption ai 7^ aj for some i, j G {0, 1, 2, 3} is rewritten that ao, ai, a2, as 
do not satisfy ao = ai = 02 = as, and it is easy to show that there exists io,ii,i2 € 
{0,1,2,3} such that ii 7^ 22, 7^ a^^ and aj^ 7^ a^j. By permutation of periods 
Ui, UJ2, UJ3 and shifts x — x + cUj (i = 1, 2, 3), we can permutate numbers ao, ai, a2, as. 
Hence it is sufficient to show the proposition under the assumption as 7^ ai and 
as 7^ a2. 

Set ci;i = 1/2, cus = r/2 and p = exp(7r-\/^r). Then ei, 62 and es are expressed as 
power series in p and we have Ci = 7r^(2/3 + 0(p^)), 62 = 7r^(— 1/3 + 8p + OijP')) and 
63 = 7r^(— 1/3 — 8p + 0(p^)). The matrix elements are expressed as 



= -4(r + 7r)(r + ^^), 4^„\,, = 64r(r + a2 - 1/2), 



(2r + a2 + as)^ - ELo ^^(^i + l)/3, 



«?r = -8{r(12r + 8ai + 12a2 + 4a3) + 47^7^ - (ai + as)^ + (ai + a2)^}. 

If p = 0, then the operator if ('o''i''2,i3) g^^^^g triangularly and eigenvalues are written 
as Ojt^ {t = 0,. . . ,d). Since the eigenvalues are quadratic in r, multiplicity of the 
eigenvalues is one or two. Hence it is sufficient to show that the eigenvalue with 
multiplicity two on the case p = separates when p varies. Assume that ai^}{= a^?|,,) 
is the eigenvalue with multiplicity two, < r < r' < d and r, r' G Z. Then we have 
r + r' = -(a2 + as). 

If r + 1 < r', then the eigenvalue around di^} is expanded as = di^} + cip + . . . , 
and ci satisfies 

(B.2) 

JfflW _f.\(7.i0) _7.{0)^(~{0) _^(0)^~(0) ~(1) 

"r,r--l"r-l,rl"r+l,r+l "r,r J J ll"r'-l,r'-l "-r' ,r' l\"'r' ,r' ''1 A"r'+l,r'+l "r',r'/ 
l"r'-l,r'-l "r',r'/"r'+l,r'"r',r'+l "r',r'-l"r'-l,r' l"r'+l,r'+l "r'.r'/i ~ 

which follows from expanding the characteristic polynomial of the matrix (af j) ^ ■_q ^ 

in variable p and observing the coefficient of p^. By a direct calculation, the condition 
that Eq. ()B.2|) for the variable Ci has multiple roots is written as (as — ai)(2r + a2 + 
as) = 0. If 2r + a2 + as = 0, then r = r' and it contradicts that r < r'. By the 
assumption as 7^ ai, it follows that Eq. ()B.2|) for the variable Ci does not have multiple 
roots, and the solution E = d^} + cip + . . . separates. 

If r + 1 = r', then r = — (a2 + as + l)/2, the eigenvalue around di^} is expanded as 
E = Or^r + Ci/2y/p + . . . , and Ci/2 is determined by 

(B.3) cl/2 = a|.*Ji,,ar^r+i = -16(ao - ai)(a2 - as)(ao + ai - l)(a2 + as - 1). 

Since < r + r' = — (a2 + as) < 2d = — (ao + ai + a2 + as), we have (1 — (ao + ai))(l — 
(a2 + as)) > 0. Combining with a2 7^ as, it follows that, if ao 7^ ai, then Eq. ()B.3|) for 
the variable Ci/2 does not have multiple roots, and the solution E = d^} + Cl/2^yp+ ■ ■ ■ 
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separates. If = oti-, then we have a^j^^ ^ = af^-^^ ^ = for r = —(0:2 + ^3 + 1)/2. The 
eigenvalue around di^} is expanded as = d^} + Cip + . . . , and Ci is determined by 

ll"r+2,r+2 "r+l,r+l A"r+l,r+l "-Ij "'r+2,r+l"'r+l,r+2 J ~ 

The condition that Eq. ()B.4j) for the variable ci has multiple roots is written as (0:2 — 
a3)(a2 + a3 — l)(2ao — 1) =0. But it is impossible, because (1 — (a2 + a3))(l — 2ao) >0 
and 03 7^ 02- Hence, if = cti^ then Eq. ()B.4|) for the variable Ci does not have 
multiple roots, and the solution E = df) + cip + . . . separates. 

Thus we obtain that the multiple roots E = df) at p = separates by expanding 
the eigenvalue as a series in p or ^/p. 

Therefore the zeros of the characteristic polynomial equation are distinct for generic 
periods (2u;i, 2ti;3). □ 

Corollary B.3. (Proposition 16'. 6|) Let /q, ^1,^2? ^3 be non-negative integers and V he 
the vector space written as Eq. (lo + h + h + h- even) or Eq. (lo + h + h + h'- 
odd). We denote the monic characteristic polynomial of the operator _f/'('o:^i.'2,^3) 
the space V by P{E). Then the roots of the equation P{E) = are distinct for generic 
periods {2lji, 2^73). 

Proof. It is shown in the proof of jHl Theorem 3.2] (or piU Proposition 3.9]) that any 
two characteristic polynomials of distinct subspaces in Eq. ()4.3|1 or Eq. ()4.8j) does not 
have common roots. Hence it is sufficient to show that the characteristic polynomial 
of the operator f/'('o:'i:'2,/3) g^^y space listed in Eq. ()4.3|1 or Eq. ()4.8|l does not have 
multiple zeros for generic periods {2u)i,2uJ3). In Proposition IB.2| it is shown that, if 
ai 7^ aj for some i,j G {0, 1,2,3}, then the characteristic polynomial does not have 
multiple zeros for generic periods. In Eqs. ()4.3l 14. 8p . the case = ai = 02 = for 
Kfo, 01,02,03 appears only for the case = ^1 = ^2 = ^3 and ao = cn = 0:2 = 0:3 = —Iq. 
We set g = lo{= h = h = h)- 

For the case Iq = h = h = h = 9, the operator H'^9,g,g,g) jg expressed as 

and the finite-dimensional space V(= y-g,-g-g,-g) for the case Iq = h = h = h = 9 
coincides with the space V for the case lo = g and = ^2 = ^3 = by replacing 
basic periods {2u}i,2uJz) — > (a;i,a;3). For the case /o 7^ and h = h = h = 0, the 
corresponding proposition is proved in Proposition IB .21 or [T^ §23.4]. Thus Corollary 
IB.3l is proved. □ 
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